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1. INTRODUCTION 

The purpose of this paper is to present a new 
method for incorporating thermal effects into 
the calculated performance of laminar lubricat- 
ing films. There is enormous interest in the 
inclusion of such effects, as the recent reviews 
of Khonsari (1,2) attest. The reason for this 
interest is well founded, since the viscosity- 
temperature dependence of typical lubricants is 
such that the viscosity can vary many fold 
across and along a bearing film, with attendant 
effects on load capacity. 

The up to date, extensive reviews in (1,2) 
remove the need here for a survey of prior lit- 
erature, and reference will be made principally 
to those works used for comparative purposes. 
Suffice to say that earlier theoretical contri- 
butions on the subject of thermohydrodynamic 
lubrication divide themselves roughly into two 
categories. In the first category are those 
which embody a full transverse (cross-film) 
treatment of the energy equation using finite- 
difference or finite-element methods, and in 
the second category are those which incorporate 
rather drastic approximations to the transverse 
phenomena, usually representing the local film 
temperature distribution by a single value. 

Both approaches certainly possess merit. But 
the first approach obtains accuracy at the 
expense of computational speed, and the second 
obtains speed at the expense of accuracy. 

We shall show here that if just two temper- 
atures, chosen at "Lobatto points'*, are used to 
characterize the transverse temperature distri- 
bution in a laminar lubricating film, the 
effects of that distribution can be surprisingly 
well predicted. The calculations we have 
so-far performed have been directed solely 
towards demonstrating this fact, and only rudi- 
mentary numerical methods have been used in the 
plane of the film. Accordingly, no meaningful 
computation times can yet be reported. But we 
believe the technique will prove to be quite 
suitable for practical calculations. 

In the present analysis, fluid properties 
are taken as constant and uniform, except for 
the viscosity. The flow is presumed to be 
laminar, with negligible inertia effects. The 
fluidity (reciprocal viscosity) is represented 
by a polynomial in terms of position across the 
film, with coefficients related to the local 
film temperature distribution. Use of this 
procedure permits a closed-form expression for 
the local lineal mass flux, albeit there is some 
difference between a fluidity profile which 
would everywhere correspond to the temperature 
profile and one which is thus approximated. 

The Lobatto-point temperatures, or mathematical 


equivalents, appear in two simultaneous partial 
differential equations obtained from the basic 
energy equation by a Galerkin procedure. 

Implementation of the present approach has 
involved considerable tedious algebra, which, 
however, once done, causes no further embarrass- 
ment. The procedure should conveniently couple 
with cavitation algorithms, and preliminary 
testing indicates that no special handling is 
required to cope with moderate recirculation at 
film entry. Moreover, it can accommodate to 
some extent the temperature streaking from hot- 
oil carryover. We therefore expect to be able 
to exploit its use in a number of interesting 
directions . 

2 . NOMENCLATURE 

Cp specific heat at constant pressure, J/kg-K 
h film thickness, m 
i Cartesian tensor index 

j Cartesian tensor index 

k thermal conductivity, J/m 2 -(K/m) 
lineal mass flux, kg/m-s 
Legendre polynomial, ith order (Pq = 1; 

Pi = C; P 2 = OC 2 - i)/2 

p pressure, Pa 

T temperature, K 

t t ime , s 

u x-wise velocity, m/s 

v y-wise velocity, m/s 

V£ i-th component of fluid velocity vector, 
m/s 

w z-wise (cross-gap) velocity, m/s 

wi Lobatto weight function for i-th quadrature 
position, Ci 

x lateral coordinate in direction of surface 
motion, m 

y lateral coordinate transverse to surface 

motion, m 

z coordinate perpendicular to gap 

midsurface, m 


Greek : 


a 

fluidity 

functions 

(see 

Table 

I) 

B 

fluidity 

functions 

(see 

Table 

I) 

4 

fluidity 

functions 

(see 

Table 

I) 

€ 

fluidity 

functions 

(see 

Table 

I) 


C dimensionless coordinate transverse to 
film, 2z/h 
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In these equations we shall treat the fluid vis- 
cosity as temperature dependent, and treat other 
fluid properties as constant. 

To these equations must be added the mass 
continuity equation for an incompressible fluid. 
Thus : 

V • v = 0 13.07} 


n fluid viscosity, N-s/m 2 

k thermal diffusivity, m 2 /s 

X fluidity functions (see Table I) 

\ fluidity (reciprocal viscosity) , m/N-s 

p fluid density, kg/m-* 

t fluidity functions (see Table 1) 

$ dissipation function, J/m 3 -s 

4> fluidity functions (see Table I) 


3. BASIC EQUATIONS 


In the absence of gravity, the momentum equation 
for a Newtonian fluid without dilational vis- 
cosity is: 



where repeated subscripts imply summation. 
And the corresponding energy equation is: 


EE + *_ 

Dt 2 \3T / Dt 
P W P 



3T 

ax. 


+ t 


13.02] 


The coordinate system used with these equat- 
ions is defined in Fig. 1. For convenience, a 
reference surface is taken midway between the 
film walls. A local coordinate system is sub- 
stituted for a fixed Cartesian system, with the 
local x-y plane tangent to this reference surface. 
The film walls are rigid, but may be moving. 

The Galerkin-style analysis used here involves 
the expansion of the temperature in a truncated 
series of Legendre polynomials. Satisfaction is 
required of as many moments of the energy equation 
as there are unknowns in this series. The ensuing 
partial differential equations for the Legendre 
components are then solved. In the present treat- 
ment, only two unknown components are used. And 
for these it is feasible to carry out explicit 
integration, as follows: 



[3.09] 


where : 


D_ 3_ 3T_ 

Dt “ 3t + V j dXj 


[3.03] 


is the time derivative following the fluid 
(Lagrangian derivative) and: 


All of the above integrals are taken from 
the "bottom** to the "top" of the film. The sub- 
scripts 2 and _2 are use< * to denote the upper 
and lower walls, respectively. The coordinate 
"z** is measured from the midplane. Continuity 
has been used to convert transverse velocity 
constructs to lateral constructs, wherever 
possible . 



is the dissipation function. 

In lubricating films, lateral diffusion of 
momentum and heat is usually much less than 
transverse. Furthermore, inertia and pressure- 
energy effects are frequently negligible, and 
the transverse variation of pressure is small. 
Therefore, we adopt the following equations for 
laminar lubricating films. 

o = - ?p + !i h !i [3 05) 

and: 


4. LOBATTO POINTS; DISTRIBUTION FORMULAS 

An expression for the lineal mass flux can be 
developed directly from Eq. [3.05]. (See, for 
example, Dowson and Hudson (3)). But the con- 
vective terms in Eq. [3.06] involve integrals 
of the velocity-temperature product, and so 
require detailed knowledge of the respective 
distributions. Such information for the veloc- 
ity is already available. For the temperature, 
we must develop our own expression. 

Consider the integral JuT dz which relates 
to the cup-mean energy flow in the x-direction. 
Let the temperatures on the two walls be known, 
and the velocity be available where needed. If 
two sampling positions for the temperature -- 
and only two — are to be allowed for estimating 
this integral, where should these positions be 
chosen? Without knowledge of end-point values. 
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it is best to choose the well-known Gaussian 
quadrature points. With knowledge of the end- 
point values, the optimum locations are the 
so-called "Lobatto points" (4), N interior 
points permitting exact numerical integration 
of a polynomial of degree (2N + 1). The case 
N = 1 yields Simpson’s rule. Here we take 
N =? 2, but it is evident that the procedure to 
be used is quite general, and that further 
refinement is possible with further analytical 
effort. 

In terms of the variable C = 2z/h, the 
Lobatto points providing exact numerical integra- 
tion over the interval -1 to +1 of polynomials 
up to the fifth degree are: 


U j T 2 + T -2 - (T 1 + T -! >} 
- T -2 - V 1 (T 1 * T -l> j I* 


[4.06] 

[4.07] 


In these expressions, the wall temperatures 
and T _2 are considered as known for purposes 
of the film calculation. It is then useful to 
note that: 


T + T 
2 -2 


[4.08] 


Location, 

Ci 

Weight, 

w i 

Subscript 

1 

1/6 

2 

1/A 

5/6 

1 

-l/A 

5/6 

-1 

l 

1/6 

-2 


Then, for example, by Lobatto numerical 
quadrature: 


/*<*« 

S' ? « 


0.4 (exact = 2/5) 


0 (exact = 0) 

0.347 (exact = 2/7) 



An expression for the fluidity might be 
developed a number of ways. Here we collate the 
fluidity with the temperatures at the Lobatto 
points; that is, at the walls and at two inte- 
rior points. The Legendre expansion for the 
fluidity is then developed in a manner com- 
pletely analogous to that for the temperature. 
For example, 

T 0 = (?_ 2 + 55_j + 5l t + 5 2 )/l2 [4.10] 

and the fluidity distribution is: 

3 

5 » £ ? v P v (f) [4.11] 

k-0 K K 


The temperature distribution Which passes 
through the Lobatto-point values is most easily 
expressed in terms of an expansion in Legendre 
polynomials. Thus, if we write: 


5. VELOCITY EXPRESSIONS; MASS FLUX 

A double integration of Eq. [3.05] (with 

S = 1/n) gives the tangential velocity vector. 

Thus: 


3 

T(C) = £ \ P k (C) 14,011 

k=0 

then the Legendre coefficients are easily 
evaluated by integration. 

£ T(C)P k CC> dc = ST*-! \ (4.021 

Or: 

iviw i4 o3) 

1=0 


v = V + A /% dC + B f a dC [5. 
-1 -1 

- - - f 1 

V 2- V -2 ' B J, « 


A = 


f 5 df 

il 


and: 


B - w 


[5.02] 


[5.03] 


The linear set of equations in [4.03] can be 
solved for the T^» providing us with two 
modes of description of the temperature distri- 
bution in the lubricating film. The T^ give 
us detail, and the T^ give us overall proper- 
ties. In particular, To is the space-mean 
temperature at the point (x,y) and *Ti is the 
first moment. 

For N 


* 2: 


(T_2 + ST.! + 5T 1 + V 

[4.04] 

- T-2 + V 1 < T 1 - T -l>j / 4 

[4.05] 


In view of the Legendre expansion for 5, [5.02] 
can be alternatively written as: 



[5.04] 


Now to obtain the lineal mass flux, the 
velocity expression [5.01] is integrated again 
across the passage, with the result: 


(V 4-V BfF* 

1 2 * -2 ; 2 3 h A 2 3 B(5 o 


[5.05] 
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Simplification gives the following more recog- 
nizable expression: 


P ■ (V 2 



[5.06] 


Here the fluidity parameters ^ are the vehi- 
cle for the ^temperature- flow interaction. The 
parameter represents asymmetry of the 

fluidity distribution. It is interesting to 
note that for symmetric temperature (and fluid- 
ity) distributions, the sole effect of temper- 
ature on the mass flow is through the arithmetic 
average of the fluidities at the Lobatto points. 
This result is a special case of the following 
formula applicable for any symmetric cross-film 
temperature distribution: 


( 5 2 + ? - 2 ) 


Vp 


ll /«■ 


5 df 


15.07] 


For such cases, Eq. [5.07] justifies the effec- 
tive viscosity concept, and shows the relative 
importance of the temperatures near the walls. 

Mass continuity applied to the mass flux 
expression [5.06] leads to a generalized 
Reynolds equation, the divergence of the mass 
flux involving spatial derivatives of pressure. 
Thus : 


£ + *.■- 0 

3t p 


[5.08] 


The first moment of the energy equation 
involves the cross velocity , w. With a little 
care, this velocity can be found via mass con- 
tinuity. We have: 

3w /3u 3 v\ ■* -» 

te-- (S**} 5 • v) z < 5 - 091 

Transforming coordinates from (x,y,z) to (x,y,C) 
we find: 


The expressions for the tangential and cross 
velocities will be essential for evaluating the 
contributions to convective heat transfer. 

6. THE TEMPERATURE EQUATIONS 

With the aid of Legendre series for the temper- 
ature and fluidity, we are now in a position to 
evaluate the integrals appearing in the zeroth 
and first moment of the energy equation; namely, 
Eqs . [3.08] and [3.09], Implementation is 
straightforward, but tedious. All requisite 
coefficients are given in Table I. 

Equation [3.08] becomes: 


-4 / T + T 

3 . :> m - 12 k 2 -2 ~ 

— (hT ) + V • “ T = ~ - T 

3t o p o h \ 2 < 


I* d£ 

PC 


- C Q [6.01] 


where "Cq" stands for the collection of 
terms: 


» ■ * £ « k T + B £ * k T' v | [6.02] 

\ k=0 k=0 

The temperature Tq appearing here is the 
C-space mean temperature, and not the mixed- 
cup temperature. 

The integral of the dissipation function is: 


/ 


2 \ ~ 2 2 ~ -» -4 

* dz - - J25 0 A + 3 ^ (2A . B) 


2 ~ 2 ~ 2 
+ I < 5 o + 5 V B 


Equation [3.09] becomes: 


a 2~ T 2 - T -2 ~ 

It <h T l ) = 60 ‘ 2 - T 1 


(6.03) 




dC - 6 Cl 


[6.04] 


where "C^” stands for the terms: 


<3 • v) z = <3 • v> c - £ 3h • 0 [5.10] 

(Here the velocity vector is understood to con- 
sist only of (u,v).) Substituting this result 
into [5.09] and integrating, we get: 


/‘I 

-> ■+ r -*■ 

V . v - f- Vh 

3v 

'h/2 | 

C h 

3C 


dz + W_ 2 

[5.11] 

Finally, integrating by parts, we obtain: 

w * ( W -2 + ? -2 * ’ 2 ) + ? 2 • - ’f * 2 * 

[5.12] 


C 1 = - hT o“-2 + ’ 


6 -2 


+ 2 ? • 1 (\? 0 + 8 £ \\ 


3 

E 

k=0 


r 


dC 


+ l E Ve • ! <K + S V 16 051 


The moment of the dissipation function is: 
f* df = | 5 X A 2 + f (5 0 + f T 2 )(2 1 • B) 


f(? 1 +f?3 ) B 2 [6.06] 
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Equations [6.01] and [6.04] are two simul- 
taneous partial differential equations in the 
two variables, and Where they appear, 

T 2 and T 3 can be eliminated via Eqs . [4.08] 
and [4.09] in favor of these dependent variables. 
Coupled with the generalized Reynolds Eq. [5.07], 
they provide our approximate thermohydrodynamic 
treatment for laminar films. 

As mentioned previously, our numerical tech- 
niques for dealing with these differential equa- 
tions have so far been extremely simple. The 
steady-state solutions to be presented in the 
next section were found by solving the foregoing 
equations in their transient form, and no study 
has been made of the possibilities for economy 
in the fluidity evaluations. Prior to this 
work, enough success has been obtained by others 
with analyses which neglect entirely any cross- 
film viscosity variations so that it seems 
unlikely that it will prove necessary to update 
the terms in Table I at every step of a solution. 

7. RESULTS FOR THE INFINITELY-WIDE SLIDER 
BEARING 

In 1963, Dowson and Hudson (3) performed some 
finite-difference calculations on the infinitely- 
wide flat slider bearing, including variable- 
property effects. They employed 100 increments 
along the length of the bearing, and a minimum 
of 20 increments transversely. These investiga- 
tors were concerned with the relative effects of 
density and viscosity variations upon load cap- 
acity, with the effects of solid-fluid thermal 
interaction, and with the possibilities for load 
support from parallel surfaces. Their findings 
serve as basis for our assumption of a constant- 
density fluid, and two special cases of their 
exploratory calculations will be used here for 
direct numerical comparisons. 

The fluid properties used by Dowson and 
Hudson were as follows (SI units): 

Density: 1.7577xl0 6 kg/m 3 

Thermal diffusivity: 7.306x10-8 m 2/s 

Viscosity: n = 0.13885*exp(-0.045(T-311,ll)) 
PaS 

Lubricant entrance temperature: 311.11 K 

Wall temperatures: Uniformly at 311.11 K 

Runner velocity: 31.946 m/s 

Bearing length: 0.18288 m 

Minimum gap: 0.00009144 m 

In the first case for comparison, the film 
thickness ratio is 2/1. Figures 2 and 3 show 
the velocity and temperature profiles obtained 
by us at the entrance, at 0.65 of the length, 
and at the exit. Figure 4 shows the correspond- 
ing pressure profile along the length of the 
bearing. Included in this last figure is the 
profile that would result if the entrance value 
of viscosity persisted throughout the film. The 
circles were read from the Dowson-Hudson curves, 
and the agreement is almost within reading 
accuracy. Figure 5 compares predicted tempera- 
ture contours. 

In the second case, the bearing surfaces 
are parallel. The possibility of load support 
through a "viscosity wedge" is being explored. 
Figure 6 compares pressure distributions. 

Again, excellent agreement is achieved. We 
note parenthetically that Dowson and Hudson 


demonstrated that, with more realistic wall 
boundary conditions for the temperature, the 
above-shown load support vanishes. 

Finally, Figs. 7 to 9 show the velocity, 
temperature and pressure profiles for a flat 
slider with 4/1 film-thickness ratio. These 
calculations were performed to test the sen- 
sitivity of the analysis to flow reversal at 
the entrance. As mentioned earlier, such 
reversal can cause problems for point-by-point 
prediction methods. We have yet to make any 
comparisons with other investigations. 

8 . CONCLUSIONS 

A Galerkin-type analysis has been made of tem- 
perature effects in laminar lubricating films. 
The procedure capitalizes on the suitability of 
so-called "Lobatto points" for sampling of the 
temperature distribution. Preliminary indica- 
tions are that the use of just two such sampling 
points enables satisfactory prediction of bear- 
ing performance even in the presence of sub- 
stantial viscosity variation. 

The procedure described herein yields two 
partial differential equations, one for the 
local space-mean temperature, and one for the 
first transverse moment of the temperature dis- 
tribution. These temperature equations are 
coupled to a generalized Reynolds equation. 
Results have been presented for the steady- 
state, inf initely-wide flat slider bearing, and 
comparisons with the earlier, detailed work of 
Dowson and Hudson are very encouraging. The 
procedure is quite general, and our intent is 
now to refine the numerical techniques, and to 
carry out calculations for more realistic con- 
figurations and boundary conditions. 
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TABLE I. - FLUIDITY FUNCTIONS 


*1 = «v 

? 2 = 5(T 2 ) 

5 -l = 

*-2 • 5 (T - 2 ) 

?0 = <* 2 + 55 x + 55 _1 + 5 _2 )/12 
h = (5 2 " 5 -2 + ^ [5 1 " 5 -l 1)/4 

h = <5 2 + 5 -2 “ 5 1 ■ 5 _1 )5/12 

? 3 “ (5 2 - «-2 1 

“0 = 2<2 ^0 - *1 + V 5) ' 3 
“l = 2( *0 " 2 V 5 + 5 3 /35)/3 

“ 2 = 2(t 0 - 2? 2 /7)/15 
*3 = 2(X 1 - 2^/37/105 

Y 0 = 2( *o - V 5)/3 

Y 1 = 2(5 ?0 " " 2? 3 /7)/15 

Y 2 = 2(25 0 - ? 2 /7)/15 
r 3 = 2(T 1 - 53 / 9 ) /35 

B 0 = 2( - 5 ?o + 4 ^1 ~ 2? 2 + 3 5 3 /7 > /15 

®j. = 2(-2t Q /5 + - T 2 /7)/3 

B = 2(f - 2? /3)/105 
.3 - 2 t 0 / 10 S 

4 0 = 2(T x /3 - 5 3 /7)/5 

4 1 = 2 <~V 5 + V 3 " *T 2 /35)/3 

4 2 = 2(4^ " T 3 )/105 

4 3 = 2 < 5 0 + t 2 /3)/35 

*0 = 0 

4> x = 2(^/3 - ? 3 /7)/5 
<t > 2 = 2(? 0 + T 2 /7)/15 
<)> 3 = 2(2T x + ? 3 /3)/I05 

£ 0 = ° 

£ i = 2( «o " V 5)/3 

e 2 = 2(^/3 - T 3 /7)/5 

L = 2? /35 
3 2 
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Figure 2.- X-Velocity versus position in gap. 
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Figure 3.- Temperature versus position in gap. 
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Figure 4.- Pressure versus position in bearing. 
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Figure 5.- Temperature contours, T/T w . 
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Figure 6.- Pressure versus position in bearing. 


TEMPERATURE, T, K VELOCITY 



POSITION IN GAP, £ 

Figure 7.- X-Velocity versus position in gap. 
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